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ABSTRACT 

We study the redshift distribution of two samples of early-type gravitational lenses, extracted from a 
larger collection of 122 systems, to constrain the cosmological constant in the ACDM model and the 
parameters of a set of alternative dark energy models (XCDM, Dvali-Gabadadze-Porrati and Ricci 
dark energy models), in a spatially flat universe. The likelihood is maximized for = 0.70 ± 0.09 
when considering the sample excluding the Sloan Lens ACS systems (known to be biased toward 
large image-separation lenses) and no-evolution, and JIa = 0.81 ± 0.05 when limiting to gravitational 
lenses with image separation A9 > 2" and no-evolution. In both cases, results accounting for galaxy 
evolution are consistent within ler. The present test supports the accelerated expansion, by excluding 
the null hypothesis (i.e., J1a = 0) at more than 4ct, regardless of the chosen sample and assumptions 
on the galaxy evolution. A comparison between competitive world models is performed by means of 
the Bayesian information criterion. This shows that the simplest cosmological constant model that 
has only one free parameter is still preferred by the available data on the redshift distribution of grav- 
itational lenses. We perform an analysis of the possible systematic effects, finding that the systematic 
errors due to sample incompleteness, galaxy evolution and model uncertainties approximately equal 
the statistical errors, with present-day data. We find that the largest sources of systemic errors are 
the dynamical normalization and the high-velocity cutoff factor, followed by the faint-end slope of the 
velocity dispersion function. 

Keywords: cosmology: cosmological parameters — gravitational lensing: strong — methods: statisti- 
cal 



1. INTRODUCTION 

In the last 15 years, several complementary ob- 
servational probes on cosmological scales have found 
strong evidence for an accelerating expansion of the uni- 
verse: distance m easurements of di st ant Type la su- 
pernovac (SNe la; iRiess et alJ (jl998| ): iPerhnutter et al.l 
(1999)), the observations of t he cosmic microwave back- 
ground anisotropies fWMAP: lBennett et a l. (2003)), and 
the baryon acoustic oscillations (BAOs) in the power 
spectrum of matter extracted from galaxy catalogs 
(jPercival et all [2001 . By assuming General Relativ- 
ity, a negative pressure component has been invoked as 
the most suitable mechanism for the observed accelera- 
tion, the simplest way of which is the cosmological con- 
stant A, with a constant equation-of-state (EoS) param- 
eter w = — 1. Dynamical dark energy models have also 
been proposed in t he lite r ature, such as the quin tessence 
dRatra fc Peebles! fl98l iCaldwell et alJ [1991 . phan- 
tom (ICaldwell et all I2002D . quintom (jFeng et al.l I2005L 
l2006t lGuo et al 1120051) . Dvali-Gabadadze - Porrati (PGP) 
(iDvali et al.l 120001: IZhu fc Alcanizl [20051: IZhu fc Serene 
2008), and Ricci da rk energy (RDE) models (|Gao et al 
2009: ILi et aHl20T(l . 



the age-redshift relation (jAlcaniz et~aT1 [20031) th e 
lookback time to galaxy clusters ( Pires et alJ 120061) . 
X-ray luminosities of galaxy cluste r s, and the Hubble 
param e ter data ([ Gaztan aga et all 120091: iStern et al 
2010t iCao. Zhu fc Lianel 120111 ICao. Liang fcZh 
201 If ). In this context, strong gravitational lens- 
ing plays an important role, providing cosmolog- 
ical tests, such as gravit a tional lensing statistic s 
(IKochanekl [l996al: [Zhu] [1991 ICoorav fc Huterer) [199 
1~ 
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The existence of a large number of theoretical models 
has triggered a variety of observational tests, based, for 
instance, on th e angular size- r edshift data of compac t 
radio sources (|Alcanizl l200"i IZhu fc Fuiimotol |2002), 
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rings in galaxy-quasar systems (|Yamamoto fc Futamasd 
I2001f) . clusters of galaxie s acting, as len s es on background 
high redshift galaxies (jSeren 3 1200a [S ereno fc Longol 
12004 ). and time delay measurements (jSchechterl 12004) . 
Results from techniques based on gravitational lensing 
are complementary to other methods and can provide 
restrictive limits on the cosmological parameters. In 
this paper we focus on one i nteresting lensing statistic 
suggested by IKochanekl (|1992D and further discus s ed an d 
developed in liter a ture (e.g., IHelbig fc Kavserl ([1996), 
IQfek. Rix fc Maozl (120031)). 

iFukugita. Futamase. fc Kasail (|1990l ) (but see also 
iNemiroffl (|1989l )') demonstrated that the gravitational 
lens expected red shift increases wi th a larger cosmologi- 
cal constant, Qa- IKocha nek (1993) derived the probabil- 
ity distribution of the redshift of the lens, as a function 
of the cosmological parameters, at fixed image separa- 
tion and source redshift. His analysis of a small sam- 
ple of four lenses strongly favored a null cosmological 
constant with respe ct to models with large J7a in a flat 
cosmological model. IHelbig fc Kavserl (|1996f) further dis- 
cussed the method, assuming no evolution of the galaxy 
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population and no detection of l ensing galaxies bey ond 
a given magnitude. Successively, iKochan ck (19965 ex ~ 
tended his analysis to include gravitational lenses with 
not measured redshift, in order to account for the se- 
lection effects. He obtained the upper limit tt\ < 0.9 
at the 95% CL, assumin g a flat cosmological model. 
lOfek. Rix k Mao3 (|2003l ) presented a new analysis of 
expected lens redshift distribution, taking into account 
number and mass evolution of the lens population, and 
applied this method to constrain both the cosmological 
and mass-evolution parameter spaces. They could ob- 
tain a strong upper limit on the cosmological constant 
(Oa < 0.89 at the 95% CL), in a flat universe with no 
lens evolution. 

Compared with other related methods based on the full 
lensing probability distribution, all of the uncertainties in 
the absolute value of the optical depth can b e elim inated 
in the differential probability of iKochaneld (|1992l ). since 
it is determined by integrating the full lensing probabil- 
ity distribution over lens redshifts. Moreover, the sharp 
cutoff of this relative probability at high redshift makes 
the quanti ty dr/dz\ (AO, a more powerful cosmologi- 
cal probe (|Ofek. Rix k Maozll200l . 

Following works investigating the evolving lens popula- 
tion have concluded that galaxy evolution is not strongly 
constrained by the redshift distribut ion test and does not 
significantly affect lensing s tatistics ([Mitchell et al1l2005l : 
ICapelo fcN ataraian 2007; Qgu ri et alj|2012j) . However, 
the evolution of mass and number density can intro- 
duce large statistical errors and bias in the analysis of 
the lens redshift distribution. Therefore, it is manda- 
tory to consider the mass and density evolution into the 
statistical analysis of the redshift distribution of grav- 
itational lenses. Other limitations include systematic 
effects due to a sample of gravitational lenses which 
completeness might be not homogenous as a function 
of the lensed image separ ation and the lens redshift 
(ICapelo k Nataraianll2007D . 

lOguri et al.l (|2012l) have presented a comprehensive 
statistical analysis of the sample of 19 lensed quasars 
found in the Sloan Digital Sky Survey (SDSS) Quasar 
Lens Survey (SQLS). This sample is used to determine 
both the cosmological constant and evolution of the mas- 
sive lensing galaxies. When considering a no-evolution 
case, a null cosmological constant is rejected at 6a level, 
providing an independent evidence for the accelerated 
expansion. 

The purpose of this paper is to extend our previous 
statistical analysis based on the angular separation dis- 
tribution of the lensed images (ICao. Zhu k Zhadl20ll 
by using the redshift distribution test to obtain novel 
constraints on the parameters of spatially flat cosmolog- 
ical models. With this aim, we use a large sample of 
122 gravitational lenses drawn from the Sloan Lens ACS 
(SLACS) Survey and other sky surveys. 

The first aim is to obtain new constraints on the cos- 
mological constant, by assessing both statistical and sys- 
tematic uncertainties, mainly due to galaxy evolution 
and sample selection. Then, we also compare a number 
of alternative dark energy models with different num- 
bers of parameter s , in our a nalys is we apply, following 
iDavis et al.l (|2007l) ; lLi etall (|2010l ). a model comparison 
statisti c, i.e., the so-c alled Bayesian information criterion 
(BIC; (|SchwarzHl97l l. 



With respect to recent works ()Oguri et al.l 120125. we 
use a larger, not homogeneous sample and focus our at- 
tention on the determination of the cosmological con- 
stant and the comparison with other alternative dark 
energy models, as the current sample size does not al- 
low a firm determination of the rate of mass and number 
evolution of massive galaxies. 

The paper is organized as follows. In Section[2j the ba- 
sics of gravitational lensing statistics is introduced, also 
allowing for number and mass evolution of the lens pop- 
ulation. We conduct a literature survey for known sys- 
tems, listing their basic parameters and defining two sta- 
tistical samples to perform the redshift test. In Section^ 
we introduce four cosmological models, and show the re- 
sults of constraining cosmological parameters using the 
Markov Chain Monte Carlo method, with and without 
galaxy evolution. In Section HJ we assess the possible 
presence of selection effects and systematic biases in our 
galaxy sample. Finally, we present the main conclusions 
and discussion in Section [5j 

2. THE REDSHIFT TEST AND THE SAMPLE 

Following lOfek. Rix k Maozl (|2003l ). the differential 
optical depth per unit redshift is 

^=n{M 1 z l ){l + z l fsJ^ 1 (1) 
dzi azi 

where n(9, zi) is the comoving number density and 
the strong lensing cross-section. 

Early-type galaxies are accurately described as singu- 
lar isothermal spheres (SIS), and it is shown that ra- 
dial mass distribution and ellipticity of the lens galaxy 
are unimportant in altering the cosm ological constraints 
(jMaoz k Rixlll99l IKochaneld Il996bl ). The SIS density 
profile is 

- 2 1 



p(r) = 



(2) 



2nG r 2 ' 

where a is the velocity dispersion and r is the projected 
distance from the center. In Section 4, we will discuss the 
systematic uncertainties introduced by this assumption. 
The corresponding strong lensing cross-section is 
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where Di, D s , and Di s are the angular diameter distances 
between the observer and the lens, the lens and the source 
and the observer and the source, respectively. Under a 
Friedman- Walker metric with null space curvature, the 
angular diameter distance reads 



D A (z ll z 2 \p) 



H (l + z 2 )J Zl E{z';p 



dz' 



(4) 



where Ho is the Hubble constant and E(z;p) is the di- 
mensionless expansion rate dependent on redshift z and 
cosmological model parameters p. The two multiple im- 
ages will form at an angular separation 



AO = 
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(5) 



We use the empirically determined velocity dispersion 
distribution function (VDF) of early-type galaxies: fol- 
lowing previous works, our sample are limited to lensing 
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early-type galaxies. The VDF i s generally modele d by a 
modified Schechter function as ([Sheth et al. 2003) 



dn 
da 



a 

— I exp 



1 



(6) 



where a, /3, n*, cr* are the faint-end slope, high- velocity 
cutoff, characteristic number density, and velocity dis- 
persion, r espectively. In the following analysis we use the 
results of lChoi et al.l (|2007|) , who analyzed data from the 
SDSS Data Release 5 to derive the VDF of early-type 
galaxies. The best-fit values of those VDF parameters 
are rt* = 8.0xl0 _3 /i 3 Mpc -3 , where h is Hq in units of 
100 km s" 1 Mpc" 1 , o-* = 161 km s -1 , a = 2.32 ± 0.10, 
and /3 = 2.67 ± 0.07. 

We also allow for evolution of the quantities n* and c* , 
by adopting the following parameterization: 

n*(z[) = n*(l + zi) p , a*(zi) = <x*(l + z\) u , (7) 

where P and U are constant quantities. 

In the following, we check the effect of redshift evo- 
lution by letting P and U be free parameters in Sec- 
tion [3j instead of adopting the evolution of (P, U) = 
(—0.23, —0- 1) pre dicted b y the semi-analytic model of 
iKang et al.l ([20051) ; iChad (|2007t) . Moreover, since the 
main goal of this paper is to constrain cosmological pa- 
rameters, we first consider P and U as free parameters, 
obtain their best-fit values and probability distribution 
function, and marginalize them to determine constraints 
on the relevant cosmological parameters of interest. 

Straightforward calculations lead to the optical depth 
per unit redshift for a system with image separation AO 
and source redshift z„. 



ia+l 



£.(A6,z s ) = T N (l + z l )l- Ua + p l 

x(l + «)»3ffl?g(#) (8) 
X exp if> (1 +*)-"/»] , 

where the normalization, 

and 

A^^Svrf^) 2 ^. (10) 
V c / D s 



For a given lens system, dr/dzi gives the relative proba- 
bility of finding the lens at different redshift: 



dzi 



dr i f Za dr ^ 
dzi J dzi 



(11) 



Besides the probability of lensing by a deflector at z\, 
dr/dzi, the total optical depth for multiple imaging r 
and the probability that a background source lensed 
by an SIS galaxy with image separation AO, dr/dAO 
(|Cao. Zhu fc Zhaol[20ll ICao fc ZhullMl) . can also be 
obtained by integrating the differential probability in 
Equation (JHJ. 



Moreover, it should be noted that the image separa- 
tion is taken into account as a prior in the method ap- 
plied in this paper, which m akes it an advan tage to use 
almost all the known lenses (jKochan ck 1992). We have 
compiled a list of 122 gravitational lenses from a vari- 
ety of sources in literature. Their basic data (lens and 
source redshifts both and the largest image separations) 
are summarized in Table [T] As mentioned above, in or- 
der to build an homogeneous galaxy sample, we limit our 
analysis to galaxies with early-type morphology. 

The main source is given by the SLACS project 
([Bolton et al.l l2008f l. providing 59 lenses in our list. 
These lenses have redshifts in the range from zi ~ 0.05 
to 0.5, making the lower redshift part of our overall sam- 
ple, with the lensed so urces ranging from z s ~ 0.2 to 
1.2 (jBolton et al.l 12008). As a consequence of the ini- 
tial spectroscopic selection method, all the SLACS grav- 
itational lenses have known spectroscopic redshifts for 
both source and lens, giving the SLACS sample an im- 
mediate scientific advantage over strong-lens candidate 
samples selected from imaging data. However, it is 
known that the SLACS sample is biased towards moder- 
ately large-separation lenses (A9 > 2"; (jArneson et al. I 
I2012D). leading to biased e stimates in the redshift test 
(jCapelo fc Nataraian|[2007l) . 

A majority of the sample was observed as part of the 
CASTLES program 0, but it also con tains gravitationa l 
lenses found in the COSMOS surve^Fa ure et al II2008T) 
and t he Extended Growth Strip (EGS; (jMoustakas et all 
12051 ). including six additional COSMOS and EGS sys- 
tems discovered recently: COSMOS5921+0638, COS- 
MOS0056+1226, COSMOS0245+1430, "Cross", and 
"Dewdrop". Finally, we also include five early- type 
gravitational lenses fr om Lenses Structure and Dy- 
namics survey (LSD. (iKoopmans fc Treul 120021 120031 : 
iTreu fc Koopmansl 120041 ) ). spanning the redshift range 
0.48 < zi < 1.00 (Q0047-2808, HST15433+5352, 
MG2016, CY2201-3201 and CFRS03.1077). 

Moreover, lenses dominated by a group or 
cluster potential will affect t he constraint re- 
sults (iKeeton fc Zabludol [2001 lOguri et al.1 [20051 : 
iFaure et al.ll2011h and previous versions of the lens red- 
shift test have certainly gone to some effort to exclude 
these systems. Therefore, systems known to be strongly 
affected by the presence of a group (e.g., B1359+154) 
or galaxy cluster (e.g., Q0957+561, SDSS1004+4112, 
and B2108+213) have been excluded from our sample. 
Moreover, we adopt the image separation criterion of 
4" to remove lenses that are influ enced by complex 
envir onments such as clusters (jOfek. Rix fc Mao9 
12003ft . On this basis we have removed Q0047-2828, 
and RXJ0921+4529 from the final sample. However, 
it should be noted that this technique is not biased by 
this criterion, because the redshift of a lens that has a 
large se paration is very low com pared with the source 
redshift llOfek. Rix fc Mao3l200l . 

The redshift distribution test requires a statistically 
complete and well characterized sample. As our list 
includes galaxies from a variety of surveys, using very 
different observational strategies and discovery spaces 
(the SDSS spectroscopic sample, Hubble Space Telescope 



4 http:/ /cfa-www. harvard.edu/castles/ 
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(HST) field sky survey, etc.), it is mandatory to verify the 
completeness of our final sample of gravitational lenses 
and its usability for the redshift test. In our analysis, we 
will use the two following samples (see Table [2]) . 

Sample A. Sixty-three lenses from the above list, ex- 
cluding the whole SLACS sample. This sample is ex- 
tracted by the same parent population as the pr imary 
sample investigated in lCapelo fe Natara ian (2003). 

Sample B. Seventy-one lenses with image separation 
larger than 2". This choice is motivated by the fact that 
the SLACS sample is biased toward moderately large- 
separation lenses and large velocity dis persions and is 
less t han 50% complete below A9 = 2" ([Arneson et al. I 

mn). 

We will discuss the samples selection functions and the 
impact on our results in Section [4[ 

3. STATISTICAL ANALYSIS 

Our statistical analysis is based on the maximum like- 
lihood technique. For a sample of Al multiply imaged 
sources, the likelihood C of the observed lens redshift 
given the statistical lensing model is defined by 

iV L 

ln£ = £ln<$p,(p), (12) 
i=i 

where Spi(p) is the particular differential probability 
given by Equation normalized to one, and p are the 
cosmological model parameters (e.g., Q\, f2 m ), for the 
Ith. multiply imaged source. Accordingly, the \ 2 is de- 
fined as follows: 

X 2 = -21n£. (13) 

and he best-fit model parameters are determined by mini- 
mizing x 2 ■ Our analys is is based on the publ icly available 
package COSMOMC (ILewis fc Bridldl2002h . 

We consider four different cosmological models to be 
tested with the observed lens redshift distribution: the 
ACDM model and three phenomenological models in 
which the vacuum energy is described as a dynami- 
cal quantity: the so-called XCDM model with the EoS 
w = pj p a free parameter, the DGP model arising from 
the brane world theory, and the RDE models. These 
models are motivated by the well known fine-tuning and 
coincidence problems of the standard ACDM model. 

We note that the previous precision cosmological ob- 
servational data have hinted that both the RDE and 
dark energy with EoS w < —1 may have du bious 
stability problems ([Feng fc Lil [20091: lAmanil 120111 ), and 
the DGP model has alrea d y been ruled out obser- 
vationally dFang et all 120081 : [D urrer fc Maartend [2010; 
iMaartens fc Koyamal 1201 Oi l, so it is indicated that these 
are just supposed to be a representative set, instead of 
viable candidates for dark energy. 

It is well known that the x 2 -statistics alone is not suf- 
ficient to provide an effective way to make a comparison 
between different models. In t his paper, we u se the BIC 
as a model selection criterion (jSchwarzl 119 78) . The BIC 
is defined by 

BIC = -2 In C max + fc In A , (14) 

where £ max , fc, N are the maximum likelihood, the num- 
ber of parameters, and the number of data points, respec- 
tively. Under this selection criterion, a positive evidence 



against the model with the higher BIC is defined by a 
difference ABIC = 2 and a strong evidence is defined by 
ABIC = 6. 

A spatially flat universe is assumed throughout the 
paper, which is strongly supported a combined 5-years 
Wilkinson Microwa ye Anisotropy Pro be (WMAP5), 
BAOs, and SN data (jHinshaw et al.ll2009[ ). As mentioned 
above, we also consider the case of an evolving popula- 
tion of lensing galaxies in order to assess the accuracy 
of our results. We consider simultaneous constraints on 
the galaxy evolution and cosmological parameters. In 
order to derive the probability distribution function for 
the cosmological parameters of interest, we marginalize 
P and U and perform fits of different cosmological sce- 
narios on both Samples A and B. Results are shown in 
Figure Q] - [6] and summarized in Table O 

3.1. The standard cosmological model (ACDM) 

In the simplest scenario, the dark energy is a cosmo- 
logical constant, A, i.e., a component with constant EoS 
w = p/ p = — 1. If spatial flatness of the Friedmann- 
Robertson- Walker metric is assumed, the Hubble param- 
eter according to the Friedmann equation is 

H 2 = H 2 [Q m (l + z) 3 + Q A ], (15) 

where Q m and £Ia represent the density parameters of 
matter (both baryonic and non-baryonic components) 
and cosmological constant, respectively. As ft — fl m + 
Q\ = 1, this model has only one independent parameter. 

We consider constraints obtained both samples defined 
above. While considering Sample A, the likelihood is 
maximized, C = £ max , for f^A = 0.70 ± 0.09 with no 
redshift evolution and fl\ = 0.73 ± 0.09 with redshift 
evolution, see Figure Q] Hereafter, uncertainties denote 
the statistical 68.3% confidence limit for one parameter, 
determined by £/£ max = exp(— 1/2). Data are consis- 
tent with the no-redshift-evolution case {P = 0, U = 0) 
at la. Specifically, the measured values of the two pa- 
rameters are P = -1.2 ± 1.4 and U = 0.22±°;2?. 

When using the Sample B, we find, in the no-evolution 
scenario, fl\ = 0.81 ± 0.05, consistent with the result 
from Sample A. When allowing for galaxy evolution, we 
find O a = 0.83 ± 0.05, P = -1.9 ± 4.6, and U = 0.16 ± 
0.30. 

In both cases, our findings are very close to the ones 
obtained from the ESSENCE s upernova survey data, 
Q A = 0.73 ± 0.04 in the flat case ([Davis et al.ll2007f ). and 
from the combined WMAP 5-year, BAO, and SN Union 
data (|Komatsu et al.l [20091 ) with the best-fit parameter 
£1a = 0.726 ± 0.015. Moreover, both samples used here 
exclude with large confidence (4cr level) the null hypoth- 
esis of a vanishing Q\, as obtained also by lOguri et al.1 
(2012) in their statistical analysis on the SQLS data, pro- 
viding independent evidence of the accelerated expan- 
sion. 

While detailed analysis on the constraints of the red- 
shift distribution test on the hierarchical models of 
galaxy evolution is beyond the scope of this work, 
we notice that these results obtained with both sam- 
ples are in broad agreement with previo us studies 
ijCapelo fc Nataraian|[2007l : lOguri et al.ll2012l ). in which 
no strong evidence for any evolution of the parameters 
U and P was found. Moreover, we find that the degen- 
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Table 1 

Sum mary of the Prope rties of the Strongly Lens ed S ystems, with the S LA CS Lenses Written in Bold . Referen ces: 1 - [Bolton et^jj ( 2008); 
2 - ITreu et al.l 00091; 3 -INewton et all i20Tlk 4 - [Ruff ct al. OoTlh: 5 -IKoopmans fe Treul 12002T); 6 - Koopmans & Trcu (2003); 7 - 

ITreu fc Koopmans! (120041); 8 - Kochanck ct al. (2000); 9 ~-lWinn et al. (2002b); 10 JGregg et ' al. (2002): 11 -Hall ct al. (2002) ; 12 
JO'Dea et all 1 119921); 13 -IWiklind fc Combes! 119951); 14 JCohen et al.l (120021)^ -ICrampton et alj 12003); 16 JFalco, Lehar. fe Shapirol 
(1997); 17 -Tonry & Kochanck (1999); 18 -Wisotzki ct al. (2002); 19 -Gregg ct al. (2000); 20 -Fassnacht & Cohen (1998); 21 
-Hagcn & Rcimers (2000); 22 -Lubin ct al. (2000); 23 -Kncib, Cohen, & Hiorth (2000); 24 -Munoz ct al. (2001); 25 -Schcchtcr ct al. 
(1998); 26 -Lehar ct al. (2000); 27 JLacv et alj ( 120021); 28 -Hcwctt ct al. (1994) ; 29 JSurdei et alj 119971) ; 30 -lEisenhardt et all 11999) ; 31 
-Lidman ct al. (2000): 32 -Tonry (1998); 33 -Myers ct al. (1999); 34 - Sicmiginowska ct al. (1998); 35 -Ratnatunga ct al. (1999): 36 
-Burud ct al. (2002b); 37 -Lehar ct al. (1993); 38 -Fassnacht ct al. (1996);" 39 -Winn ct al. (2002a); 40 -Morgan ct al. (2001); 41 
JLangston eTaTT (U 98911; 4 2 -Wiklind & Combes ( 1996); 43 -IWinn et all (12000 ); 44 -ISvkes et all (119981) ; 45 -Fassnacht ct al. (1991); 46 
feurud et all 120023); 47 - IHuchra et al . (1985): 48 -llnada et all (120051); 49 -IBig gs et all 120031): 50 - IJohnston ~ al. (2003); 51 
-Blackbur ne et all 120081): 52 -ISluse et all MXm ; 53 -lEigenbrod eta l. (2006); 54llnada ct al. ( 2008]) ; 5 5~IMorgan et alj 120041): 56 
-IQfek et al (2006); 57 -IFaure et all 120081); 58 {Anguita ct al. (2009); 59 -Lagattuta ct al. (2010); 60 -Moustakas ct al. (2007). 
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Table 2 

Summary of the Subsamples Used in the Analysis of This Paper. 



Subsample 
Sample A 
Sample B 
Sample C 
SLACS 

Full sample 




P 

Figure 1. Simultaneous constraints on the cosmological constant 
and VDF evolutions in the flat ACDM model obtained by using 
Sample A. Upper panel: likelihood distributions as a function of 
Qa with and without redshift evolution. The red dotted line is ob- 
tained after marginalizing over P and U. Lower panel: constraints 
on P and U after marginalizing over f^A- Dashed lines in the lower 
panel represent the case without redshift evolution (P = U = 0). 




°Jio-8 -6 -4 -2 2 4 6 8 10 

P 

Figure 2. As in Figure [l] using Sample B (i.e., lensing systems 
with image separation larger than 2"). 

erate direction in the evolution parameters corresponds 
to a constant lensing probability, as already been noted 
bv lQguri et all (I2012L 

Previous studies on lensing statistics (|Chael I2003t 
lOfek. Rix &: Maol 12003ft considering a not-evolving ve- 
locity dispersion function have obtain ed consistent re - 
sults with the galaxy number counts dim et al.l 12002ft . 
iMitchell et ail (|2005f ); iCapelo k, Nataraianl (|2007T ) as- 
sumed a non-evolving shape for the VDF and obtained 
results consistent with earlier results. To sum up, all 
previous results on redshift evolution from strong lens- 
ing statistics appear to be in agreement with no-evolution 
of early- type galaxies. As we only considered early- type 
lensing galaxies, our results further confirm this conclu- 
sion. 

Substantial evolution of galaxy number (i.e, P ~ — 1) 



Definition 



= 0.8 

63 lenses from Table [1] excluding the SLACS sample 

71 lenses with image separation larger than 2" Q ~J 

51 lenses with image separation not larger than 2" ^* 

59 lenses from the whole SLACS sample — — Q.6 

122 lenses from Table[T] !Q 
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and mass (i.e., U ~ 0.25) are suppo rted by studies of 
early- type galaxie s from z = 1 to (jFaber et al.l 120071 : 
iBrown et all 120071; lOguri et all 120121) This conclusion is 
consistent with our results in Section [3] 

3.2. Constraints on selected dark energy models 

Now we focus on the constraints obtained on selected 
dark energy models. Here, we only use Sample A, as 
results obtained with Sample B are coherent at la level. 

3.2.1. Dark energy with constant equation of state (XCDM) 

While deviating from the simple case w = —1, the 
EoS of dark energy w < —1/3 can also give birth to 
an accelerated university expansion. In a zero-curvature 
universe, the Hubble parameter reads: 

h 2 = Hl\n m (i + zf + n x (i + z) 3 ^} . (i6) 

Obviously, when flatness is assumed (f2 TO + Q\ = 1), it 
is a two-parameter cosmological model, p = {Cl x , w}. 

The best-fit values of the parameters are: Sl x = 
0.77 ± 0.17; to = -2.3tl$ witn n0 redshift evolution 
and tt x = 0.79 ± 0.13; w = -2.1^1 with redshift evolu- 
tion, see Figure [3] and Figure H] for the confidence limits 
in the Q x — w plane. However, we note that the lower 
limits on the parameter w are probably an artifact of 
the prior w > — 5, which may be tested and constrained 
with future larger lens sample. Also in this cosmologi- 
cal scenario, the lens redshift data are consistent with no 
redshift evolution: when marginalizing over fl x , we find 
P = -1.4 ± 1.4 and U = 0.20 ± 0.28. The Einstein's 
cosmological constant (to = — 1) is still consistent within 
la. Meanwhile, compared to the cosmological constant 
model, this flat cosmology with constant EoS dark energy 
provides a lower Xmim but a higher BIC: ABIC = 1.72 
with no redshift evolution and ABIC = 1.55 with redshift 
evolution. However, we note that comparing with the 
cosmological constant model, the two-parameter XCDM 
model performs relatively well under the information cri- 
terion test. 

3.2.2. Dvali-Gabadadze-Porrati model (DGP) 

In a brane world theory, the accelerated expansion of 
the universe can be explained by the gravity leaking 
out into the bulk at large scales with a crossover scale 
(jDvali et al.l l2~000). The Friedmann equation is 

3M|j (V - f") = + zf , (17) 

where Mp\ is the Planck mass and r c = (H (l — f2 m )) _1 
is the crossover scale at which the induced 4-dimensional 
Ricci scalar dominates. In this model, the Hubble pa- 
rameter expresses as 

H 2 = flg(v/fi m (l + zf + Ct rc + v^T) 2 , (18) 

and the DGP model has one one free parameter p = 
{f2 m } in a flat universe. 

For the DGP model, the best-fit parameters are Vl m = 
0.25l°;Jg with no redshift evolution and ft m = 0.22t°;Jg 
with redshift evolution (See Figure [S]). We find that the 
DGP model is somehow worse than the ACDM model 
with this observational data. While its Xmin * s l ar g er 








Figure 3. Likelihood contours for the flat XCDM model at 68.3% 
and 95.4% CL in the Q x — w plane obtained by using Sample A 
with no redshift evolution. The lower limit on w is probably an 
artifact of the prior w > — 5 
and the horizontal dotted line indicates a cosmological constant 
with w = — 1. 

than that of the ACDM model by about 1.2, it gives 
ABIC = 1.20 with no redshift evolution and ABIC = 
1.19 with redshift evolution. 

3.2.3. Ricci dark energy (RDE) model 

Following [Gao et al.l (|2009l ): lLi etall (|2010l) . the aver- 
age radius of the Ricci scalar curvature |7£| -1 / 2 might 
provide an infrared cutoff length scale. Accordingly, the 
DE energy density is 

p de = Zp 2 {H + 2H 2 ), (19) 
and the Hubble parameter reads 

H 2 = ^[^(l + -) 3 + (1 - fr|)(l + (20) 

where /3 is a positive constant to be determined. This is 
a two-parameter model with p = {r2 m , /?}. 

For the RDE model, the best-fit values of the parame- 
ters are fi m = 0.22^°;^ ; P = 0.29±0.19 with no redshift 

evolution, and Q, m = 0.18tjH^ /3 = 0.28 ±0.18 with red- 
shift evolution. In Figure [6l we plot the likelihood con- 
tours in the Q m — f3 plane. Given the largest information 
criterion result: ABIC = 1.78 with no redshift evolu- 
tion, and ABIC = 1.62 with redshift evolution, the RDE 
model performs the worst in all the cosmological models 
considered in this paper. 

4. BIASES AND POSSIBLE SYSTEMATIC EFFECTS 

Thus far we have considered only statistical errors. In- 
deed, cosmological tests based o n strong lensing h ave 
been somehow controversial since iKochanekl (|1996a[) . m 
particular for the possible biases ass ociated with sample 
selection ([Capelo fc Nataral an 2007). In this section, we 
discuss several possible sources of systematic errors, in- 
cluding: sample incompleteness, unknown survey selec- 
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Figure 4. As in Figurc[T] but for the fiat XCDM model obtained 
by using Sample A with redshift evolution. The lower limit on w 
is probably an artifact of the prior w > — 5 
and the horizontal dotted line indicates a cosmological constant 
with w = — 1. 

tion function, uncertainties in the lensing galaxy proper- 
ties and lens modeling, in order to verify their effect on 
the cosmological constraints. 

First of all, one general concern is given by the fact that 
strong gravitational lenses are a biased sample of galax- 
ies. Most of the previous works found no evidence for a 
biased sample of the lensing popula tion with respect to 
massive early- type g alaxies (see, e .g., Jjreu et all l|200fih ). 
On the other hand, iFaure et al.l (|201ll ) found possible 
evidence for the stellar mass of lensing early-type galax- 
ies to evolve significantly with redshift. However, it is 
still not clear whether this supports a stronger lensing 
bias toward massive objects at high redshift or if it is a 
consequence of the possible higher proportion of massive 
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a 

m 

Figure 5. Normalized likelihood as a function of Q m for the 
DGP model obtained by using Sample A with and without redshift 
evolution. 

and high stellar density galaxies at high redshift. This 
could be addressed in dedicated numerical simu lations 
(|van de Ven et alJl2009tlMandelbaum et al.H2009h . as the 
available lens samples cannot allow yet to discriminate 
between the two alternatives. 

We now estimate the systematic errors due to statisti- 
cal sample incompleteness. As both our Samples (A, B) 
have been put together from different surveys, differences 
in the observing strategies (and selection functions) may 
cause systematical errors hard to estimate. In order to 
evaluate the effects due to a selection bias, we have red- 
erived the best estimate on tt\ (with a no-evolving lens 
population) by using the whole catalog of n — 122 lenses 
and two additional, smaller samples (see Table [5]): the 
first one includes only galaxies from the SLACS, and the 
second one the gravitational lenses with separations not 
larger than 2", including 59 and 51 systems, respectively. 
All of these three samples clearly suffer from strong se- 
lection effects or are very inhomogeneous and are there- 
fore not suitable for deriving constraints by means of the 
redshift distribution test. However, they can shed light 
on the amplitude of the possible systematic errors due 
to sample inhomogeneity, incompleteness, and selection 
bias toward non-representative systems. For instance, it 
is well known that the SLACS catalog is characterized by 
a selectio n function favoring mo derately massive ellipti- 
cals (e.g., lArneson et al. I (|2012l )): also the somehow ex- 
treme case of Sample C is not complete (below A6 = 2") 
as the probability to detect a lensing galaxy is related 
to the image separation as very small-separation lenses 
easily escape detection (both in present-day imaging and 
spectroscopic surveys). Therefore, we expect that the es- 
timates of fl\ obtained from these samples allow us to 
establish an upper limit on the systematic errors due to 
a not well-defined lens sample. 

Results are shown in Figure [7] When using the appar- 
ently incomplete sample of small separation lenses (Sam- 
ple C) we find n A = 0.65lg'i5. This rather smaller value 
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Figure 6. Likelihood contours for the RDE model at 68.3% and 
95.4% CL obtained by using Sample A with no redshift evolution 
(upper) and with redshift evolution (lower). 

can be related to the fact that lower-mass lenses, pro- 
ducing small-separation images, will tend to be located 
at redshifts lower than expected in a large-f^A model. 
Hence, this determines a slightly lower value for the cos- 
mological constant. 

When adopting the whole SLACS sample, we find 
£!a = 0.73^0 is- Finally, by considering the whole, inho- 
mogeneous list of 122 lenses, we obtain fl A = 0.71+ u ;g£. 
When we compare these values with the best estimate 
obtained from Sample A, we note that systematic errors 
do not exceed ^0.1 on the cosmological constant. 

We have assumed all the lenses to be isolated sys- 
tems, with negligible line-of-sight contamination. How- 
ever, current studies fi nd that proximate galaxies 
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( Cohn fc Kochanekl I2004D and environmental groups 
Keeton fc Zabludofil I2004I) can have various effects on 



"a 

Figure 7. Normalized likelihood plot as a function of Oa for the 
flat ACDM model, for the three biased samples described in the 
text. 

the primary lens galaxies, even though the most signif- 
icant effects only appear to be biasing the galaxy cllip- 
ticities and image multiplicities. Nevertheless, most of 
lenses in Sample B come from the SLACS su rvey where 
the ro le of environment has been assessed in iTreu et all 
(12009? ). Namely, it was found that for SLACS lenses the 
typical contribution from external mass distribution is 
no more than a few percent. Therefore, the environmen- 
tal effects on observed image separations appear to be 
relatively small (certainly smaller than statistical errors 
arising from the current sample size of lenses) . 

In addition to the main lens galaxy, the contribution 
from line-of-sight density fluctuations to the lens poten- 
tial should also be taken into c onsideration . Bas ed on 
the final lens sample from SQLS. lQguri et aLl ()2012[ ) have 
investigated the line-of-sight effect and found that its ef- 
fect on the total lensing probability is rather small com- 
pared with the contribution of other system atic errors to 
the sy stematic error on (see Table 2 in lOguri et al.l 
(2012)). Meanwhile, it is noted that this effect may 
also have a larger impact on lensing probabili ties for the 
lens samples with larger images separations (lOguri et al.l 

Evolution of the VDF and model uncertainties can in- 
troduce additional systematic errors. Here we estimate 
these systematic errors on the constraint results of the 
flat ACDM with the f ull sample (n = 12 2 lenses), in a 
similar way as done in lOguri et al.1 ((2012). The analysis 
in Section [3] suggests that unconstrained redshift evolu- 
tion of the velocity function is one of the most significant 
sources of systematic error. An additional source of un- 
certainty is the relation between velocity dispersions and 
image separations. This uncertainty is not only related 
to the difference between the velocity dispersion asis 
of the mass distri bution and the observ ed stellar veloc- 
ity dispersion <7n (| White fc Davislll996f) . but also many 
complexities such as the velocity dis persion normaliza - 
tion factor for non-spherical galaxies (jOguri et al.ll2012l) 
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Figure 8. Upper panel: constraint results on the cosmological 
constant and the parameter for the flat ACDM model with 
the full sample (n = 122 lenses). The horizontal dashed line indi- 
cates the fiducial value Je = 1 assumed in the paper. Lower panel: 
constraints on with different values of f& considering a fiducial 
error of ~ 5% on the values of image separations. 

and the detailed luminosity profiles of galaxies. 

Hence, we introduce the parameter Je that combines 
the velocity dispersion asis and the stellar velocity dis- 
persion <7o 

cs/s = /Wo, (21) 

The Ei nstein radius given by Equation [TO] is then modi- 
fied as (IKochaneklll99l IQfek. Rix fc MaxTz|[2003l) : 

^-■"(t)*!^- (22) 

More specifically, we have kept Je as a free param- 
eter to mimic the following effects of: (1) system- 
atic errors in the difference between ctq and ers/s; 
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Figure 9. Normalized likelihood plot for the flat ACDM model 
with the full sample (n = 122 lenses), by introducing the uncer- 
tainties on the VDF parameters: the faint-end slope a and the 
high- velocity cut-off /3, as listed in Tabic [3] 



2) the influence of line-of- sight mass contamination 
Keeton. Kochanek. fc S cliak 1993); (3) softened isother- 
mal sphere model which may change the typical image 
separations (jNaravan fc Bartelmannlll996| k and (4) the 
effect of secondary lenses on observed image separation. 
Recently, by combining stellar kinematics (central ve- 
locity dispersion measurements) with lensing geometry 
(Einstein radius determination from position of images) , 
iCao et all <|2012f ) have tested a combined gravitational 
lens data set including 70 data points from SLACS and 
LSD, and obtain ed constraints on fp cons istent with the 
previous results (IQfek. Rix fc M aoz 200l . 

We consider the flat ACDM model and obtain simul- 
taneous constraints on £7a and Je with the full sample 
(n = 122 lenses). The constraint result is shown in Fig- 
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Table 3 

Summary of the Priors and Standard Systematic Allowances 
Included in the Analysis. 



Table 4 

Constraint Results Obtained by the Full Sample (n = 122 lenses) 
for the Flat ACDM Model with Different Systematic Errors. 



Parameter 


Allowance 


EoS of DE (w) 
DE density in XCDM (il x ) 
DE density in ACDM (fi A ) 
Matter density (f2 m ) 
Evolution of n* (P) 
Evolution of cr„ (U) 
Normalization factor (f E ) 

Faint-end slope (a) 
High-velocity cut-off (/3) 


-5 < w < 
< Q x < 1 
< n A < 1 
o < n m < i 

-10 < P < 10 
-1< U< 1 
(0.5) 1 / 2 < f E < (2.0) 1 / 2 
a = 2.32 ±0.10 
P = 2.67 ±0.07 



Systematics 



ureEJ with the marginalized constraints Je — 1.08±0.14 
and n A = 0.65±g;||. It is obvious that the cosmological 
constraint on J7a becomes much weaker, though the con- 
cordance ACDM is still favored within la error region. 
A vanishing cosmological constant is still ruled at la. 

The slightly larger value of Je indicates that larger 
fl\ is preferred by the lens statistics with /e = 1 used in 
this pap er. This tend e ncy is consistent with the previous 
result of lQguri etaLI (|2012f) and may help explain why 
the observed imag e separation seems to be h igher than 
model predictions (jCapelo fc Nataraianl2007T) . However, 
we find that the = 1 case still consists with the data 
better to than 2a. To be more specific, in order to take 
into account the measurement uncertainty and the ap- 
proximations we are doing (no lens ellipticity accounted 
and no external shear), we estimate a fidu cial error of 
~ 5% on the values of image separations (jGrillo et al.l 
2008), which is equivalent to a ~ 5% uncertainty on the 
parameter /J,. 

The velocity distribution function given by Equation [6] 
is another important source of systematic error on the 
final results. While adopting the best-fit values of the 
VDF measureme nt in the SDSS Data Release 5 by 
IChoi et all ((20071) as our fiducial model, we investigate 
how the cosmological results are altered by introducing 
the uncertainties on a and j3 as listed in Table [3l With 
the results of measurements on the VDF, we vary the pa- 
rameter of interest while fixing the other parameters at 
their best-fit values. For example, based on the n = 122 
sample, we vary the faint-end slope a by ±0.10 and find 
that this effect is quite negli gible when compared to the 
present accuracy of the test (Mitchell et al. 20051). 

The complete set of standard priors and allowances in- 
cluded in the analysis of the above systematics is summa- 
rized in Table [3] As shown in Figure [5] and and Table 0] 
by comparing their contribution to the systematic error 
on J7a for the flat ACDM with the full sample, we find 
that the largest sources of systemic error are the dynam- 
ical normalization Je and the high- velocity cutoff /3, fol- 
lowed by the faint-end slope a of the VDF. This fi n ding i s 
consistent with the earlier results in IQguri et al.l (|2012f ). 
Meanwhile, this result is also compatible to the concor- 
dance cosmological model (Oa ~ 0.7, and Vt m ~ 0.3). 
Indeed, current samples of lenses do not allow to dis- 
criminate between an J7a and a dynamical dark energy 
component. 

5. CONCLUSION AND DISCUSSION 
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JE 
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J E 


= 1.00; a = 


2.32 


P 


= 2.74 


n A = 


= 0.76 ±0.07 



Q A 



Since the discovery of the accelerating expansion of 
the universe, in addition to the standard ACDM cos- 
mological model, a large number of theoretical scenarios 
have been proposed for the accele ration mechanism. Ex- 
amples include the q uintessence dRatra fc Peebles! 1 19881 : 
ICaldwell etaD 11998 ), phantom dCaldwell et al.1 12002) , 
auintom dFeng et al.lf2005l 120061 iGuo et al.|2005l) and the 
Chaplygin gas dKamenshchik et al.l I2001L iBentq et al.1 
Alam et al. ll2003HZhull2004UZhang fc Zhul l2006) . 



Unap l 

Hill 



All these acceleration mechanisms should be tested with 
various astronomical o bservations s uch as SNela, WMAP 
(|Komatsu et al.l I2009D and BAO (jPercival et al.ll2007[ ). 
However, it is still important to use as many as differ- 
ent observational probes to set constraints on the cosmo- 
logical parameters. In this work, we have followed this 
direction and used the redshift distribution of two well- 
defined samples of lensing, elliptical galaxies drawn from 
a large catalog of 122 gravitational lenses from a variety 
of surveys (see Table [Tj). The BIC is also applied in this 
analysis to assess various dark energy models and make 
a comparison. 

Considering the two cases with and without the red- 
shift evolution of the velocity function of galaxies, we 
have analyzed four dark energy models (the ACDM, the 
XCDM with constant w, the DGP and the RDE models) 
under a flat universe assumption. For each model, we 
have calculated the best-fit values of its parameters and 
found its ABIC value. Results are plotted in Figure [TJ51 

The fit and information criteria results are summarized 
in TableO It is shown that for the zero-curvature ACDM 
model, the likelihood is maximized, for fl\ = 0.70 ± 
0.09 with no redshift evolution and fl\ = 0.73 ± 0.09 
with redshift evolution, when using the Sample A (see 
Section 2). Consistent results (within la) are derived 
when the alternative Sample B is considered. 

We have also derived simultaneous constraints on the 
redshift evolution of the parameters n* and cr» of the ve- 
locity function. The constraints in the P — U plane also 
indicate that the data are consistent with the no-redshift- 
evolution case (P = 0; U = 0) at la, with P = -1.2±1.4 
and U = 0.22±gJ|. Both in the no-evolution and galaxy 
evolution scenario, we rule out with high confidence a 
vanishing cosmological constant, with both the lens sam- 
ple s (at confidence lar ger than ~ 4<r), as recently found 
bv IQguri et al.1 (|2012l ) by adopting a smaller and inde- 
pendent lens sample. Therefore, the redshift test adds 
an independent evidence to the accelerating expansion. 

The obtained likelihood distributions shown in Fig- 
ure [T] are also in agreement with the results from an- 
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LCDM XCDM DGP RDE 



Figure 10. Graphical representation of the values of ABIC 
for each model, relative to the LCDM (ACDM). The horizontal 
green line indicates ABIC=2, which is considered positive evidence 
against the model with the higher BIC. 



Table 5 

Summary of the constraint results and information criterion with 
Sample A. The ABIC values for all other models are measured 
with respect to the ACDM model. 



Model 



Constraint result 



ACDM (P = U = 0) 


n A 


= 0.70 ±0.09 









n A 


= 0.73 ±0.09 


_ 9 q+1-3 
— z.0_ 2 7 





XCDM (P = U = 0) 




= 0.77 ±0.17; w 


1.72 


(P + + 0) 


Q x 


= 0.79 ±0.13; w 




1.55 


DGP (P = U = 0) 




- 25+ 011 




1.20 






= 0.22±°;J° 




1.19 


RDE (P = U = 0) 




= 0.22™ ,3 = 


0.29 ±0.19 


1.78 




f2 m 


= 0.18±° : ^;/? = 


0.28 ±0.18 


1.62 



ABIC 



alyzing data of WMAP 5-year results with the BAO 
and SNe Union data, and the larg e-scale structures in 
the SDSS luminous red galaxies (iSpergel et all 120031 : 
iTegmark etaD l2004t lEisenstein et al.l 120051 ). which im- 
plies that gravitational lensing statistics provides an in- 
dependent and complementary support on the ACDM 
model. 

We give a graphical representation of the BIC test Fig- 
ure [10] Following the ACDM model, the DGP model is 
the only one-parameter model that reduces to the ACDM 
but gives a worse fit, as we obtained close values for the 
matter density. The XCDM model gives a comparably 
good fit, but has one additional free parameter, that is 
in accordance with the best-fit ACDM model (within la 
range for the EoS parameter). The other two-parameter 
model, RDE, provides a worst fit to the data, though 
the difference in BIC (ABIC) indicates no clear evidence 
against it. Therefore, while still not firmly ruling out 
competing world models, the redshift distribution test 
clearly favors the cosmological constant m odel, a con- 
clusio n in accordance with previous works ([Davis et al.l 
[2001 . 



In order to asses the reliability of our results and the 
related systematic errors, we have introduced three ad- 
ditional lens samples, characterized by different degrees 
of inhomogeneity and incompleteness, and re-derived an 
estimate for the cosmological constant assuming a not 
evolving lens population. Two lens samples (the whole 
heterogeneous catalog of 122 systems and the SLACS 
sample) give results fully consistent with those discussed 
above, confirming us that systematic errors due to sam- 
ple selection are not larger than statistical uncertainties. 

Our model involves several uncertainties and assump- 
tions that introduce additional systematic errors in our 
cosmological analysis (see Table [3]). By comparing the 
contribution of each of these systematic errors to the sys- 
tematic error on the flat ACDM model (see Figure [8] and 
[9] and Table [4} , we find that the largest sources of system- 
atic error are the dynamical normalization and the high- 
velocity cutoff factor, followed by the faint-end slope of 
the velocity dispersio n function, which i s consistent with 
the earlier results in lOguri et al.l ([2012D . Moreover, the 
comparable systematic errors suggest the importance of 
careful studies of the systematics for robust cosmological 
constraints from lens statistics. 

Finally, we note that four important effects should 
be mentioned. Evolution of the source population can 
also matter the technique applied in this paper, which 
might have a sma ll second-order effect on the statistics 
([Qguri et al.1 120121 ). We have also neglected systematic 
uncertainties due to the effect of small-scale inhomo- 
geneities on large-scale observations. In fact, the inho- 
mogeneous matter distribution can affects l ight propaga- 
tion and th e related cosmological distances ([Sereno et al.l 
120011 120021 ) . a lthough its effect on the total lensing statis- 
tics is small (jCovone et al.l [20051 ) . Meanwhile, the lens 
redshift test applied in this paper may also lose the sta- 
tistical power of absolute lensing probabilities. The last 
one is, multiple errors or biases in the method could eas- 
ily be canceled out, which may lead the result to be a 
statistical fluke (|Maozll2005l ). 

Despite some of its inherent difficulties, the red- 
shift distribution test, with either larger gravitational 
lensing samples from future wide-field surveys such 
as Pan-STARRS and Large Synoptic Survey Tele- 
scope by taking advanta ge of time-domain information 
(|Qguri fc Marshalll l2010h or a joint investigation with 
other cosmological observations, could be helpful for ad- 
vancing such applications and provide more stringent 
constraints on the cosmological parameters. 
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